Creating stable molecular condensate using a generalized Raman adiabatic passage 

scheme 
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We study the Feshbach resonance assisted stimulated adiabatic passage of an effective coupling 
field for creating stable molecules from atomic Bose condensate. By exploring the properties of the 
coherent population trapping state, we show that, contrary to the previous belief, mean-field shifts 
need not to limit the conversion efficiency as long as one chooses an adiabatic passage route that 
compensates the collision mean-field phase shifts and avoids the dynamical unstable regime. 
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Molecules offer a whole new dimension in the study 
of ultracold atomic physics. In particular, resonant 
photo- or magneto- association (Feshbach resonance pj) 
of cold atoms to molecules represents a matter- wave ana- 
log of second harmonic generation and has become a new 
paradigm of coupled macroscopic quantum systems. Due 
to energy conservation, such processes generally produce 
molecules in a vibrationally and/or electronically excited 
quasibound level and hence are not energetically stable. 
This does not seem to be a serious problem in the case 
of fermionic atoms due to the suppression of molecular 
decay by Pauli blocking. In fact, several groups have now 
successfully Bose condensed so-formed molecules Q. In 
contrast, although evidence of macroscopic coherence has 
been observed in several experiments, excited molecules 
formed by a pair of bosonic atoms have very short lifetime 
(< a few ms) and the atom-molecule conversion efficiency 
is limited to < 10% jUQ. Hence it is very important to 
be able to create deeply-bound ground state molecules 
from atomic Bose condensates with high efficiency, which 
is the focus of the current work. 

First proposed for photoassociating nondegenerate 
atoms into stable molecules 0, stimulated Raman adi- 
abatic passage (STIRAP) aided by Feshbach resonance 
is considered a more efficient way of converting atomic 
condensates into molecular ones than the bare STI- 
RAP implementation of photoassociation 0, @ • In this 
scheme, the free atomic, the quasibound and the ground 
molecular states form the three-level A-system to which 
STIRAP can apply 10]. The success of STIRAP relies on 
the existence of the coherent population trapping (CPT) 
state . In a linear A-system, the CPT state exists when 
the two-photon resonance condition is satisfied, hence 
STIRAP can be straightforwardly implemented by ap- 
propriately choosing the laser frequencies. In the case of 
condensate, collisions between particles give rise to non- 
linear mean-field shifts which dynamically changes when 
population is transferred from atomic state to molecular 
state. This poses a serious problem for STIRAP as col- 
lisions shift the system out of the two-photon resonance. 
For typical experimental parameters, a conversion effi- 



ciency of only ~ 20% is predicted with a somewhat com- 
plicated laser sequence containing seven Raman pulses 
[5j . A possible remedy is to use a low-density atomic con- 
densate where the effect of collisions can be minimized 
[TT|. This, however, creates new problems as stronger 
coupling fields and longer time scale are required. 

The purpose of our work is to show that, despite of the 
presence of the mean-field collisions, the "two-photon" or 
the CPT condition [See Eq. @ bellow] can be dynami- 
cally maintained and high atom-molecule conversion rate 
can be achieved in condensates of typical density. The 
nonlinear collisions, however, may induce dynamically 
unstable regimes in the parameter space. It is crucial for 
the success of STIRAP to avoid these unstable regimes 
when designing the route of adiabatic passage. 




FIG. 1: The energy diagram of three-level atom-molecule sys- 
tem involving free-quasibound-bound transitions. 

Our model system is schematically sketched in Fig. ^ 
where the relevant energy levels are denoted by \a) (free 
atom), \m) (quasibound molecule) and \g) (ground state 
molecule). Levels \a) and |ra) are coupled by a magnetic 
field through Feshbach resonance with coupling strength 
a' and detuning e (e is experimentally tunable via ex- 
ternal magnetic field), while \m) and \g) are coupled by 
a laser field with Rabi frequency and detuning Q and 
A, respectively Without loss of generality, we take 
both a' and fl to be real as their phase factors can be 
absorbed by a trivial global gauge transformation of the 
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field operators. In the interaction picture, the Hamilto- 
nian describing the system reads 

H = hJ dv\ ^A^J(r)#t( r )^.( r )^( r ) 



+ e*L(r)*m(r) + — |*]„(r)tf (r)tf a (r) + h.c. 



(A + e)M(r)* fl (r) 



n r 



*lW* s (r) + /i.e. 



where "I^ (i — a, m and g) is the annihilation field op- 
erator for state the terms proportional to Ay- rep- 
resent two-body collisions with X' u = Airhai/rrii and 



Ay = X'ji = 2tt ha,ij I for i ^ j (a^ and ay are s-wave 
scattering lengths, and /Uy is the reduced mass between 
states i and j) characterizing the intra- and inter-state 
interaction strengths, respectively. Here we consider a 
uniform system and hence have dropped the usual ki- 
netic and trapping terms. 

From the Hamiltonian we can easily derive the equa- 
tions of motion of the field operators. We adopt the 
standard mean-field treatment to replace the field oper- 
ators ^i/^/n with c- number tpi, where n is the density 
of the total particle number. The mean-field approach 
ignores high order quantum correlations, but is expected 
to be valid for systems with a sufficiently large number 
of particles |l3| . The set of the mean-field equations are 



dtpg 

dt 
dt 
dt 



~i (A a \lpa\ 2 + Km \lpm\ 2 + X ag \lpg\ 2 ^j "Pa ~ iaip m 1p*, 

-11pm ~ ieipm - i (A m \lp m \ 2 + X am \lp a \ 2 + X mg l^gl 2 ) 1p m - + i^lpg, 

-1 (Xg \lpg\ 2 + X a g \ ifj a \ 2 + X m g | ?/> m | 2 ^ ij) g ~ 1 (A + e) l\) g + £ — ^Wl, 



(la) 
(lb) 
(1c) 



where Xi = X' u n, Ay = Ay-n, a = ol \fri are the renor- 
malized quantities, and the term proportional to 7 in 
Eq. i|lb(l is introduced phenomenologically to simulate 
the decay of the quasibound molecules. Implied in 
Eqs. is the assumption that this decay dominates all 
the other particle loss mechanisms. Hence, Eqs. Q hold 
only in a time scale much shorter than the characteristic 
times associated with other loss mechanisms, although 
the time scale itself can be much longer than 7 . 

The technique of STIRAP allows the system to evolve 
in a coherent superposition of stable particle states (|o) 
and \g)) over time, effectively eliminating the loss of par- 
ticles through highly unstable states (|m)). It is impossi- 
ble to apply the STIRAP in a usual sense directly to our 
model because our model contains nonlinear interaction 
terms. Another difference with the conventional STI- 
RAP is that while the coupling strength of the |m) <-> \g) 
transition can be controlled by the optical pulses, the 
coupling strength of the \a) «-> |m) transition a is fixed 
not by optical means but by the hyperfine interaction 
and is hence independent of time [l4| . 

We now proceed to show that despite these important 
differences, Eqs. Q support a CPT steady state with 
ip m = 0. To this end, let us first neglect the particle 
loss by taking 7 = and treat all the other parameters 
as time-independent. Then the total particle number is 
conserved, i.e., \ip a \ 2 + 2 \i/j m \ 2 + 2 \ip g \ 2 — 1, so that we 
can introduce the atomic chemical potential /i through 
the following steady-state ansatz 

ipa = \ip a \ e M -e-^*, ip m , g = l^mj e^e'^. (2) 



Putting into and taking the time derivative to be 
zero, one finds that the following CPT solution exists: 

|^| 2 = 2 =1-2 1^| 2 , (3a) 

yd + 8 (a/O) 2 + 1 

K| 2 = 0, (3b) 

with 9g = 26 a and the corresponding chemical potential 
/i = A a j ip® j 2 + X ag j tpg j 2 , under the condition 

A = -e + (2A ag - Xg) K| 2 + (2A a - A ag ) |^| 2 . (4) 

Equation 0} is the generalized "two-photon" reso- 
nance condition which incorporates the nonlinear colli- 
sional phase shifts. This CPT solution is thus a general- 
ization of the one found in the collisionless limit y| ■ 

Two crucial conclusions regarding the utility of this 
CPT state in the adiabatic atom-molecule conversion can 
be arrived by observing (J3J and (@J : i) Population is con- 
centrated in states |a) and \g) under the respective limit 
a/Q — > and a/VL — > 00, which facilitates adiabatic co- 
herent population transfer between the atoms and the 
ground state molecules by tuning the optical Rabi fre- 
quency f2 (a, as we mentioned, is regarded as fixed in 
time); ii) Since the temporal dependence of O uniquely 
determines the time evolution of the population, we can 
design a priori the temporal evolution of the laser de- 
tuning A in accordance of £l(t) such that the resonance 
condition Q is maintained dynamically at any time. (A 
similar idea can be found in a study of vortex coupler 
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|l5| . ) In so doing, we eliminate the need to minimize 
the effect of collisions by reducing the density of the sys- 
tem ^lj or by other means |Tfjj . 

The existence of the CPT state, however, does not 
guarantee that this state can be followed adiabatically. 
We have to study the stability properties of the state. 
This has been neglected in all the previous studies on this 
problem. For this purpose, we adopt the linear pertur- 
bation approach, adding small fluctuations (in both am- 
plitudes and phases) to the steady-state CPT solutions 
and linearizing the equations of motion. The excitation 
frequencies of the system other than the zero frequency 
Goldstone mode can be found analytically as 
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u = ±J[b± \/b 2 -4c) /2, 



(5) 



where 

b = n 2 /2 + 2a 2 \i/j° a \ 2 +A\ 

c = (tt 2 /4 + a 2 |^| 2 ) 2 + y |^| 4 ( A 9 ~ + 4A ") A ' 

with A = A | 1 2 + ^ag - 0.5A mg - e and A = 2A Q - \ ag - 
A Q m + 0.5A mff . When u> becomes complex, the corre- 
sponding CPT state is dynamically unstable. Hence the 
unstable regime is given by either c < or c > 6 2 /4. The 
instability here is caused by collisions as it can be easily 
seen that in the absence of collisions, lo is always real. 
The typical results from the stability analysis based on 
the parameters of our interest are summarized in Fig. 
where the e)-space is divided into the stable (white) 
and the unstable (dark) regions. 

In our calculations, we have taken the parameters for 
23 Na. The s-wave scattering length for sodium atom is 
well-known and is around a a = 3.4 nm |17| which yields 
\' aa = 1.18 x 1CT 16 m 3 s" 1 . We choose a' =4.22 x 10" 6 
m 3/2 g -i^ correS p 0nc Ji n g to the atom-molecule coupling 
strength for the sodium Feshbach resonance at a mag- 
netic field strength of 85.3 mT jig. We take the con- 
densate density n to be 5 x 10 20 m~ 3 . This gives rise to 
a = 9.436 x 10 4 s" 1 , and A a = 5.9 x 10 4 s" 1 = 0.625a. 
There are so far no good estimates on molecular scatter- 
ing lengths. So we take the collisional coefficients involv- 
ing the molecular levels to have the same magnitudes 
as A a . Generally speaking, there are two unstable re- 
gions: Region I is thin along e-dimension and centers at 
e = A a g — X mg /2 in the small limit and at e = 2A — A am 
in the large f) limit; Region II occurs at small SI and lies 
either above or below Region I in e-direction depending 
on whether A g — 4A ag + 4A Q is positive or negative. We 
will show that it is very important to avoid these unstable 
regions in order to convert atoms into stable molecules. 

Samples of our results on atom-molecule conversion are 
displayed in Fig- EI where we have numerically solved the 
full set of Eqs. including the loss term with a time- 
varying Rabi frequency given by 




n(t) 



1 — tanh 



t-tr 



(6) 



1 2 3 4 5 6 



FIG. 2: Stability diagram in (fi, e)-space. S7 and e are in 
units of a. Other parameters: X a = 0.625, A m = A 9 = 0.1875, 
X am = X ag — \ mg — 0.1875, all in units of a. 



and A changes accordingly as determined by JIJ. Al- 
though in principle the Feshbach detuning e can be var- 
ied temporally to optimize the conversion efficiency, we 
find it not quite necessary. Fixing e at a finite value not 
only simplifies the experimental procedure, more impor- 
tantly it also allows us to avoid strong condensate losses 
near the exact Feshbach resonance (e = 0) ^ n ex_ 

perimental implementation of our scheme may thus take 
the following steps. First, turn on the coupling field and 
fix its amplitude at Jl max while the Feshbach transition is 
far-off resonance. Next, we bring the Feshbach transition 
suddenly into near resonance. This is followed by the adi- 
abatic ramp of the amplitude and the frequency of the 
coupling field according to Eq. © and Eq. respec- 
tively. This implementation represents a counterintuitive 
sequence as the optical field coupling the initially empty 
state \g) and \m) is turned on before the interaction on 
the transition involving the initially populated state \a). 

Figure Ofa) shows the time evolution of the popula- 
tion \4>a,g,m(t)\ 2 ■ Here e is chosen such that the system 
remains in the stable regime (see Fig. Also plotted 
in the figure is the analytical CPT solutions of Eqs. [3] 
As it can be seen, the exact population dynamics follows 
closely the prediction of the CPT solutions. The small 
discrepancies at later times can be attributed to the fact 
that the system can not maintain adiabaticity completely 
[20|. The final population in the stable molecular state 
is \ip g (oo)\ 2 — 0.415, corresponding to a conversion ef- 
ficiency of 77 = 2\ipg(oo)\ 2 — 83%. The conversion ef- 
ficiency can be further improved if the optical field is 
ramped down more slowly. In contrast, in the dynamics 
depicted in Fig.^b), initially the populations follow the 
CPT solutions, but significant deviation starts to occur 
at about t = 200/a when the system enters the unstable 
regime. The final conversion efficiency is less than 20%, 
which cannot be improved by using a more slowly varying 
optical field as one does in the stable regime. Figure[Sfc) 
summarizes the conversion efficiencies for various e. The 
marked asymmetry between positive and negative e is a 
dramatic manifestation of the dynamical instability. 
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FIG. 3: (a) and (b): population as functions of time for e = 
—3 and 3, respectively. Also plot in (a) are the CPT solutions 
of Eq. ®. Here 7 = 1, fW = 40, r = 40, to = 120, and 
all the other parameters are the same as in Fig. |2] Time is 
in units of 1/a, all other quantities are in units of a. (c) 
Conversion efficiency r\ — 2\ip g (t = oo)| 2 as a function of e. 

In conclusion, we have presented here an effi- 
cient method to convert atomic condensate into stable 
molecules using a single optical pulse without the need 



to sweep through the Feshbach resonance. Our method 
generalizes the conventional STIRAP scheme in that the 
laser is frequency chirped to compensate for the mean- 
field shift arising from the particle collisions, so that, at 
least in principle, the "two-photon" resonance is main- 
tained at all time. We also point out that the presence of 
particle collisions also gives rise to potential dynamically 
unstable regimes in parameter space. A high conversion 
efficiency can be reached only if one designs the STIRAP 
route such that these unstable regimes are avoided. 

Note that we do not expect the presence of the trap 
will change our results significantly so long as the trap- 
ping frequencies are made sufficiently small. An estimate 
in the Thomas- Fermi limit and for a spherical trap of fre- 
quency ujq — 2ir x 41 Hz indicates that more than 90% 
of a total of 5x 10 6 sodium atoms will occupy the central 
part of the trap where the variation of the mean-held 
shift due to density inhomogeneity is less than a/10. A 
numerical calculation based on the same parameters as 
in Fig.[3]but with the two-photon resonance condition vi- 
olated by an amount of a/10 indicates that a conversion 
efficiency as high as 76 % is still attainable. 

Finally we want to mention that although we specif- 
ically studied a Feshbach resonance situation, due to 
the formal resemblance of the governing equations, our 
method can be straightforwardly applied to a photoasso- 
ciation situation. In fact, we solved Eqs. (3) of Ref. [HI 
with our laser frequency modulation scheme and ob- 
tained a final stable molecular population of 0.42 at den- 
sity n = 4.3 x 10 14 cm~ 3 using the parameters for their 
Fig. 2(c). Compared to 0.08 obtained in Ref. [ll|, the 
advantage of our method is quite obvious. 
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